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ABSTRACT 

We construct a real space renormalization group (RG) approach for Ising spin glasses on 
hypercubic lattices within the scheme of the Migdal-Kadanoff approximation using replicas. 
Our replica symmetric solution yields results consistent with simple decimation previously 
obtained and the introduction of breaking of replica symmetry within the RG is discussed, 
which inserts in a natural fashion non-linear RG into the problem. 



PACS. 75. ION - Spin glass models 
PACS. 75.50L - Spin glasses 
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1. Introduction 

The study of the magnetic properties of the so-called spin glasses (random magnetic 
systems with quenched disorder and competing interactions) has been challenging the imag- 
ination of many workers in the field, and up to now the only fairly complete theory is the one 
associated with the solution of the mean field Sherrington-Kirkpatrick (SK) model (for an 
excellent account of the development of the field see Fisher and Hertz 1991), which suffers 
from the original sin of having size dependent and infinite range interactions. Nevertheless, 
these limitations of the model might be forgiven when one considers bona fide models and 
from them obtain the known mean field theories either for the pure ferromagnetic case (Kac 
1968) or spin glasses (de Almeida 1992). 

Usually the mean field results hold at very large dimensionalities, give the qualitative 
aspects that a theory in finite dimensions should encompass, and are obtained without much 
effort. Spin glasses changed this altogether for even its mean field theory (assumed as the 
solution of the SK model) took many years to be worked out, and only in the past few years 
a phase transition in Ising 3d systems has been largely accepted (sec Fisher and Hertz 1991). 

To study the nature of the phase transitions in 3d systems a natural framework would 
be to employ the renormahzation group techniques (see Wallace and Zia 1978 for a review). 
However, from the beginning the absence of a correct picture of the nature of the cooperative 
fluctuations in spin glasses will be hanging over any results which might be obtained. 

Among the various forms of implementing the renormahzation group approach, real space 
methods in lattice models has been widely used for its conceptual and technical simplicity, 
in particular the Migdal-Kadanoff renormahzation group (MKRG) approximation. In this 
work we study Ising spin glasses on hypercubic lattices using the MKRG approximation 
in conjuction with the replica method. The study of Ising spin glasses using the MKRG 
approximation is not new (Ney-Nifle and Hilholst 1993 and references therein). What is 
new in our work, we believe, is the introduction of the replica method and new parameters 
allowing breaking of replica symmetry. Thus a full consistence between the mean field (SK) 
theory and a renormahzation group approach may emerge. 
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2. The Model and its RG 



Let us consider a model described by the foUowing Hamihonian 



7^ = — ^ JijaiUj 
{hi) 



±1 



(1) 



where labels the pair of nearest neighbor sites on a hypercubic lattice with random 
exchange interactions coupling constants Jj^ 's assumed to have zero mean and variance J of 
a gaussian distribution. The bond average of the replicated partition function of the model 
(1) is (Edwards and Anderson 1975) 



where < ••• > means bond average, K = (3 J, and a, (3 = l,2,...n are replica indices. 
Let us now perform an approximate RG transformation on (2) within the Migdal-Kadanoff 
approach. Here we carry this out by working with the effective Hamiltonian in replica space 
as defined by (2): 



As far as we know, none of the existing works on spin glasses using the RG approach used 
(3) as the starting point. First, consider the simple decimation scheme using a rescaling 
factor b — 2, and a new rescaled lattice and interactions are obtained from the previous ones 
using the clusters indicated in figure 1. Calling Si and Sj the spins of the rescaled lattice 
each connection in the cluster shown in figure 1 will produces a contribution to the effective 
interaction between Si and Sj, which is obtained by tracing out the replicated intermediary 
{a°') spins in the following expression 






(3) 





Using standard transformation equation (4) assumes the form 





where 



3 



Ki(xi,yi) = tgh-^ [tgh(Kxi)tgh{Kyi)] 



(6) 



and 



(7) 

From (5) it is seen that each connection in figure 1 contributes with a term Ki{xi^yi) to the 
efi^ective interaction between Sf and 5'°', all replicas with the same weight, each term with 
a probability distribution specified by (7). The new effective interaction is then given by 



2d-i 



(8) 



1=1 



for d-dimensional hypercubic lattices. The probability distribution of the interactions in 
the decimated lattice is then given by 



P{K') = / , 



K'-J2 tgh-' {tgh{Kx,)tgh{Ky,)) 

1=1 

whose mean deviation < K'"^ > is 

<K'^> ^1-JdzJ dK' jJ{Dx^Dy^{K'fexp{izK'- 
- iz^tgh'^ [tgh{kxi)tgh{Kye)] \ 
e ) 

which may be recast into the form {h = 2) 



(9) 



(10) 



d-l 



Dx Dy \tgh~^ {tgh{Kx)tgh{Ky)) 



(11) 



Analogously, one obtains for the fourth moment of the distribution 



< K'^ >= 3 (< K'^ >y + 2'^-^ 



< li* > -3 < >2 



(12) 



where < > is the double integral in equation (11) and 

<u'^>= j Dx j Dy [tgh-^ {tgh{Kx)tgh{Ky))f (13) 

Equations (11) and (12) could have been obtained without the need of introducing the 
replicas shown in equation (4). We exhibit them to stress its replica symmetric nature. 
They are the result of a one step decimation starting from a gaussian bond distribution 



and the last term in (12) is a measure of the deviation from a gaussian of the new iterated 
distribution. They may be extended to provide an iterative map between the moments of the 
effective bond distributions at each RG transformation. Here we take the approximation 
of considering (11) and (12) vaUd at each RG step which amounts to approximating the 
iterated probabihty distributions by gaussians, which is a good approximation as d assume 
large values. 

Prom (11) we obtain that there is a phase transition only for d > 2, and Tc ~ 0.95, 2.12, 
3.51, 5.32 for d = 3,4,5, and 6, respectively. The correlation coefficient p =< K'^ > /3(< 
K''^ >Y, which is one for a gaussian distribution, takes the values at T^, p = 1.11 {d = 3); 
1.11 {d — 4); 1.08 {d — 5) 1.05 {d — 6) justifying the above mentioned approximation. The 
correlation-length exponent (i/) defined by the flow away from the non-trivial flxed point of 
(11) is given by 

and we find u = 1.41 {d = 3), u = 0.77 {d = 4), z/ = 0.60 {d = 5), u = 0.55 {d = 6). For 
comparison the high temperature series expansion for these exponents (Klein et al 1991) are 
u = 1.37 (d = 3), 1/ = 0.95 (d = 4), i/ = 0.73 {d = 5) for a ±J distribution. Others critical 
indices may be obtained in the usual way invoking their relationship obtained through usual 
scahng. 

It might be worth mentioning that phenomenological scaling (Bray and Moore 1985, Bray 
1988 and references therein) claim that the ordered state is governed by a zero-temperature 
fixed point, such that the width of the renormalized probability distribution at low temper- 
atures has a growth characterized by an exponent y, i.e., K' ~ 2^^. Prom (11) we find 
numerically that |/~((i — 5/2)/2, yielding a lower critical dimension dc~ 2.5. For compar- 
ison, numerical studies (Bray and Moore 1985, Bray 1988) estimate y ^ —0.30 {d = 2) and 
~ 0.2 {d — 3) while our result gives y ~ —0.25 {d — 2) and y ~ 0.25 {d — 3). However, we 
note that our results thus far have been obtained within the context of replica symmetric 
approach. An improved treatment may upset the above fiow of the renormalized probability 
distribution. While equation (11) gives a stable zero-temperature fixed point for d > d^, a. 
proper account of frustration may turn it unstable as in highly frustrated systems studied 
by McKay et al (1982). 



We consider now the question of introducing breaking of replica symmetry within the 
present approach. The renormahzation group framework for studing critical properties can 
rarely be carried out exactly. In general, one has to resort to some approximations and 
frequently introduce some parameters to optimize the results (sec Wallace and Zia 1978). 
In the case under study, we may take this step writing equation (4) in the form 



where the gala's are new parameters introduced into the theory to be determined say, for 
instance, variationally from a free energy. Breaking of replica symmetry (RSB) within the 
real space RG may then be pursued by iterating (15). A first step RSB scheme (see below) 
leads to a iterative map involving several parameters which is under study. 

In this work we have only analysed the replica symmetric case = q = 1 and leave for 
future work other cases. 

3. Discussion 

In this work we have reported, we believe, the outline of the renormahzation group in 
real space as should be applied to spin glasses (here only Ising systems were considered but 
extension to other systems is straightforward). Only the naive replica symmetric case was 
considered and results consistent with previous works were obtained (Southern and Young 
1977, Curado and Meunier 1988, Bray 1988, Ney-Nifle and Hilhorst 1993, and references 
therein). It was our aim in this work to explicitly point out the replica symmetric (RS) nature 
of the results, obtain only the simplest results within the RS solution, and to lay the ground 
for extension of the theory. We are not aware of any similar discussion except in a work of 
one of the present authors (de Almeida 1993), where it was argued that it is in the nature of 
these complex systems (spin glasses) that they require the use of nonlinear renormahzation 
group in their analysis. This may completely change standard universality rules as obtained 
for uniform systems (critical exponents depending only on the space dimension and on the 
number of spin components) which have been confirmed numerically through Monte Carlo 
calculation (Bernardi and Campbell 1993) for o? = 3 and 4, using several distinct forms 




(15) 
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for the probability distribution of the interactions. Moreover, a study of the multifractal 
spectra of the order parameter for short range Ising spin glasses in d — 3 (Coutinho et al 
1993) indicate a nontrivial /(a)-function for T <Tc thus suggesting the complex nature of 
a proper scaling for spin glasses. In the present approach, nonlinear renormalization group 
sets in imediately as soon as one start breaking replica symmetry among the Qa/s in equation 
(15), in a fashion similar to the breakdown of linear response theory and the Fisher relation 
in the spin glass phase (Bray and Moore 1980). 

In the last few years many works either experimental (Lefioch et al 1993 and references 
therein) or numerical (Badoni et al 1993, Hetzel et al 1993, and references therein) have 
argued in favour of the mean field SK picture in finite dimensions. There is even the sug- 
gestion that for Ising spin glasses there is no qualitative difference in their behaviour for 
all dimensions above the lower critical one (Bhatt et al 1991). Our preliminar analysis of 
breaking replica symmetry using equation (15) following a Blandin's like seminal scheme 
(Blandin 1978), indicates y exponents much larger than in the RS case suggesting the pos- 
sibility of existence of an de Almeida-Thouless transition hne (de Almeida and Thouless 
1978) within the present context. In addition, RSB introduces an unfolding of the bonds 
probability distribution much like in line with the RSB picture yielding a probability law 
for the probability law of the overlap of the order parameters (Mczard et al 1984). All these 
topics deserve further investigation which we leave for future work. In this letter we hope 
to have added another step forward in the long and hard spin glass way. 
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Figure Captions 

Figure 1: Migdal-Kadanoff b — 2 cluster. There are p — 2*^"^ connections joining the S^s 
spins in d-dimensions. 
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